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ABSTRACT 

We use the clustering of galaxies around distant active-galactic nuclei (AGN) to derive an estimate of 
the relationship between galaxy and black hole mass that obtained during the ancient quasar epoch, at 
redshifts 2 ;$ z ;$ 3, when giant black holes accreted much of their mass. Neither the mean relationship 
nor its scatter differs significantly from what is observed in the local universe, at least over the ranges 
of apparent magnitude (16 ^ Gab ~ 26) and black-hole mass (10 6 ;$ Mbh ~ 1O 1O5 M0) that we are 
able to probe. 

Subject headings: galaxies: high-redshift — cosmology: large-scale structure of the universe — quasars: 
general 
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1. INTRODUCTION AND RESULTS 

The study of black holes has been driven to the fore- 
front of extragalactic research by the recent discovery 
of black holes as massive as a billion suns inside nearby 
bulge galaxies. Simple physical arguments (e.g., Silk & 
Rees 1998) suggest that these enormous objects should 
profoundly affect the process of galaxy formation, a be- 
lief that is strengthened by the tight observed correla- 
tion between the masses of local galaxies and their black 
holes (Gebhardt et al. 2000; Ferrarese & Merritt 2000). 
Various theoretical models attempt to explain the exis- 
tence of the correlation with a wide range of physical 
processes. Since these models make discordant predic- 
tions for the evolution of the correlation over time, we 
decided to test them by measuring the relationship be- 
tween galaxy and black hole mass in the distant past, at 
redshifts 2 < z < 3. 

A novel approach (see, e.g., Kauffmann & Haehnelt 
2002) let us use our existing surveys (Steidel et al. 2003; 
Steidel et al. 2004) of ~ 1600 galaxies at redshifts 
1.5 £j z ;$ 3 to measure the dependence of galaxy mass 
Mh on black hole mass Mbh over a 5-decade baseline 
of black hole mass, reaching masses roughly 1000 times 
smaller than the limits of other surveys (e.g., Shields et 
al. 2003; Walter et al. 2004; Croom et al. 2005) at sim- 
ilar redshifts. After using the technique of Vestergaard 
(2002) to estimate the masses of the black holes that 
powered each of the 79 active-galactic nuclei (AGN) in 
our survey (see figure and the appendix), we estimated 
the typical halo mass for black holes in different mass 
ranges by measuring how strongly the other galaxies in 
our survey clustered around them. 

Adelberger & Steidel (2005) describe our analysis in 

1 Based, in part, on data obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership between the Cal- 
ifornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
of the W.M. Keck Foundation. 
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Fig. 1. — Overview of the characteristics of the AGN in our 
sample. Upper left: Redshifts and absolute AB magnitude at rest- 
frame 1350A. The uncertainty in the AB magnitude is < 0.2 mag- 
nitudes for even our faintest objects (e.g., Steidel et al. 2003). 
Upper right: Relationship between CIV line width and apparent 
AB magnitude at rest-frame 1350A. The uncertainty in line width 
ranges from 10—20%, and is dominated by systematics (e.g., con- 
tinuum placement) for the brightest AGN. Lower panels: Relation- 
ship between CIV line width, mi3so , and the resulting estimate of 
black- hole mass Mbh- The selection bias is severe in our AGN 
sample, since (for example) we deliberately targeted AGN that 
were bright and had broad emission lines. These panels show the 
characteristics of our sample as selected, not of a fair sample of 
high-redshift AGN. 



more detail. Briefly, we estimated the cross-correlation 
length ?*o from the number of galaxy-AGN pairs with 
angular separation 60" < 8 < 300" and comoving ra- 
dial separation AZ < 30/i _1 Mpc with the approach of 
Adelberger (2005), then used the GIF-LCDM numerical 
simulation (Kauffmann et al. 1999) of structure forma- 
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Fig. 2. — Relationship between cross-correlation length ro and 
implied mass of the black-holes' galaxies at redshifts 2 ^ z 3. 
Blue, green, and red points show the dependence of ro on Mh 
in the GIF-LCDM simulation (Kauffmann et al. 1999) at redshifts 
z = 2.97, 2.74, 2.32. Small offsets have been added to the abscissae 
for clarity. The relationship is somewhat uncertain because it de- 
pends on the assumed masses of the halos that contain non-active 
galaxies. Error bars show the la uncertainty in the relationship 
that results from the uncertainty in the masses of the non-active 
galaxies. 



tion in a standard cosmological model to estimate from 
tq the total (dark matter plus baryon) mean mass 
of the galaxies associated with black holes in each mass 
range. The relationship between ro and Mh depends on 
the redshift and on the mass of the typical (non-active) 
galaxies in our survey, but the resulting systematic errors 
in Mh are small compared to the random errors (figure^ ■ 

We found galaxy-AGN cross-correlation lengths of 
r = 5.27t\H for the 38 AGN with 10 5 8 < M BH /M Q < 
10 s and r = 5.20±i;f| for the 41 with 10 s < M W /M Q < 
lO 10 5 . The inferred relationship between log(MBH) and 
(log(M/j)) is shown in Figure|3| 

If the predicted relationship between galaxy and black- 
hole mass has the form log Mbh = /(logMft) + e, with 
/ a function to be specified and e a random deviate, 
the expectation value of log(Mh) for a given value of 
log(MBH) follows from the elementary relationship 
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(1) 



where l h = lag(M h ), l B n = log(M BH ), P(k) is the dis- 
tribution of log(Mh) measured in the GIF-LCDM sim- 
ulation and extrapolated with the appropriate Press- 
Schechter (1974) formula, and P(Ibh\Iii) is the distribu- 
tion of log(MBH) at fixed galaxy mass, which depends 
on / and on the characteristics of the random variable 
e. Solving equation ^ numerically for different functions 
/ under the assumption that e has a Normal distribu- 
tion with rms a e , we find the theoretical tracks shown in 
Figure 01 

The solid blue line is for a M BB \-Mh relationship iden- 
tical to the one observed locally, log(MBH/10 7 M Q ) = 
1.651og(AVlO 12 Af ) + e (Ferrarese 2002). For this line 
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Fig. 3. — Observed and expected relationship between black hole 
and halo mass (in solar units) at redshift z ~ 2.5. The ordinate is 
(log(.M h )| log(AfBH)): the mean value of log(Af h ) for a given value 
of log(MBir). Points show our observations. Vertical error bars 
show the ler random uncertainty. Horizontal error bars show the 
mean and rms value of log(AfaH) for the two groups of black holes. 
Lines show theoretical predictions. Solid blue line: no evolution 
in mean M^n-M^ relationship, negligible intrinsic scatter. We 
adopt equation 6 of Ferrarese (2002) for the local relationship, but 
her two alternatives fit our observations comparably well. Short- 
dashed blue line: no evolution, 1 dex of scatter in Mbh at fixed 
Mh- Green line: local relationship scaled by (1 + z) 5 / 2 . Red 
line: relationship at z = 3 in a supply-limited accretion model (Di 
Matteo et al. 2003). Observations at large A/bh agree with any 
of these scenarios, as has also been noted by Shields et al. (2003) 
and Walter et al. (2004). Our small- A/bh data help distinguish 
between them. 



we assumed tr £ = 0.5, roughly the expected error in our 
black hole masses (Vcstcrgaard 2002). The line therefore 
assumes negligible intrinsic scatter in the correlation. It 
fits the data well. 

The other lines show that alternative relationships in 
the literature generally provide a worse fit. The green 
line results from scaling the ratio of black hole to galaxy 
mass by (1 + z) 5 / 2 , as advocated by many semi-analytic 
models (e.g., Haehnelt, Natarajan, & Rees 1998; Wyithe 
& Loeb 2002; Volonteri, Haardt, & Madau 2003). The 
red line shows the redshift z — 3 prediction A/bh/Mq = 
6.2 x IO^M^/IO^Mq) 1 033 of a model in which black 
holes accrete a fixed fraction of the total gas mass in 
each merger (Di Matteo et al. 2003). The dashed blue 
line assumes that the mean Mbh-M;, relationship is the 
same as observed locally but that its intrinsic scatter has 
increased to 1.0 dex. Increasing the scatter decreases 
the typical mass of galaxies that contain black holes of 
a given mass. This is because galaxies with low masses 
are much more common than galaxies with high masses; 
when the scatter in the Mbh-M;, relationship is big, 
the largest black holes are more likely to reside in low 
mass galaxies with unusual ratios of Mbh/M^ than in 
high mass galaxies with normal ratios. The clustering of 
galaxies around AGN would therefore be far weaker than 
we observe if there were no relationship at all between 
Mbh and Mh- 

A x 2 test suggests that the three alternatives to the no- 
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evolution model ((1 + z) 5 / 2 scaling, supply limited accre- 
tion, large a e ) disagree with the observations at the 90- 
95% level. They can therefore be considered marginally 
consistent with our present data, although the odds are 
against them and more extreme evolution from the local 
relationship (e.g., Haehnelt & Rees 1993) can be ruled 
out with high significance. 

The apparent lack of evolution in the Mbh-M^ corre- 
lation seems consistent with models in which the corre- 
lation results from active feedback from the black hole. 
In these models the black hole mass is pinned near the 
maximum allowed by its halo at all times. If this max- 
imum is set by the escape velocity at a fixed proper ra- 
dius from the black hole, it will not depend strongly on 
redshift. One might object that black holes are able to 
enter the quasar phase in these models only because their 
masses have temporarily fallen below the maximum al- 
lowed by their growing halos, and so the most luminous 
AGN should never lie on the correlation. As long as the 
quasar phase occurs near the end of the accretion, how- 
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ever (e.g., Hopkins et al. 2005), the black hole should 
have nearly achieved its equilibrium mass. In any case, 
a slight decrease in Mbh at fixed Mh would make the 
predictions fit our data even better. 
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APPENDIX 



TECHNICAL DETAILS 

The general population of galaxies tends to cluster more strongly around individual galaxies with larger masses. We 
exploit this effect to estimate the masses of the galaxies that harbor black holes. After estimating the characteristic 
mass M g of the general galaxy population from its measured correlation length (Adelberger et al. 2005), we use the 
GIF-LCDM simulation to calculate as a function of M/, how strongly galaxies of mass M > M g cluster around galaxies 
of mass M > Mh- We infer the masses of the galaxies that harbor various black holes by finding the value of Mh 
required to match the observed cross-correlation length ro- Figure [5] shows the relationship we used to estimate from 
our measured cross-correlation length ro the typical mass of the galaxies containing the black holes (green points). 
Adopting other plausible relationships between ro and Mh would change the inferred masses by less than their random 
uncertainties. Percival et al. (2003) and Kauffmann & Haehnelt (2002) have shown that halos undergoing mergers 
have the same correlation length on large scales as other halos of the same mass, so our estimates of ro should provide 
reasonable estimates of the halos masses even if AGN are fueled by mergers. 

To estimate the random uncertainty in ro , we took a Monte-Carlo approach that exploited the similarity of the AGN- 
galaxy cross-correlation length to the galaxy-galaxy correlation length. We generated many alternate realizations of 
our data by treating randomly chosen galaxies in each field as that field's AGN, rather than the true AGN themselves, 
and recalculated ro for each simulated sample. Since the galaxies in our survey outnumber the AGN by more than 
twenty to one, the simulated samples are nearly independent of each other and of the true sample. We took the rms 
spread in ro among them as the la uncertainty in our measured measured correlation length rQ bs . The distribution of 

X 2 = 2Zt( r o bs — r o IOd )/ cr n>] 2 f° r tne predicted values of ro in figure should be roughly equal to the distribution of x 2 
in the simulated samples around the line rQ red = constant = rp S , where r$ s is the galaxy-galaxy correlation length in 
our sample. We used this distribution to associate our measured values of x 2 with a P-value. 

Our conclusion depends on the assumption that the estimated black hole masses Mbh are not wildly inaccurate. 
We estimate Mbh from an AGN's luminosity I = XL\ at A = 1350A and CIV line-width FWHM with the relationship 
that is observed in the local universe: M BH /M Q ~ 10 6 - 2 (//10 44 ergs- 1 ) - 7 (FWHM/1000kms -1 ) 2 (Vestergaard 2002). 
Correcting for a stellar contribution to the AGNs' luminosities (which we have not done) would decrease our lowest 
observed values Mbh even further, strengthening our conclusions. Our estimated black-hole masses would be too low 
for some AGN with small Mbh if their observed CIV emission line were produced in the narrow-line region rather 
than the broad- line region (as we assume). In this case the line widths would be roughly equal to the galaxies' stellar 
velocity dispersions (Nelson 2000), at least for radio-quiet AGN, but in fact the galaxies' mean stellar velocity width 
(~ 200 km s _1 ) is an order of magnitude smaller than the mean AGN line-width for Mbh < 10 8 M Q (2100 km s _1 ) 
or Mbh > 10 8 M Q (4900 km s _1 ). It is far smaller than even the smallest observed AGN line- width in our sample, 
800 km s _1 . Radio- loud AGN make up too small a fraction of our sample to affect our results if omitted. In any case, 
the observed range of Mbh is so large that our estimates of Mbh would have to be wrong by ~ 1 order of magnitude 
to alter our results significantly. We cannot rule out the idea that the relationship between A/bh, luminosity, and 
line- width was utterly different in the past, but it seems easier to believe that the relationship between Mbh and Mh 
has not changed at all. 
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